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ABSTRACT 

The Spitzer Survey of Stellar Structure in Galaxies (S G) is an Exploration Science Legacy Program 
approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d<40 
Mpc, |b| > 30°, rciBcorr <15.5, D25 >1') survey of 2,331 galaxies using the Infrared Array Camera 
(IRAC) at 3.6 and 4.5/mi. Each galaxy is observed for 240 s and mapped to > I.5XD25. The 
final mosaicked images have a typical la rms noise level of 0.0072 and 0.0093 MJy sr -1 at 3.6 and 
4.5/im, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at 
H3.6tim{AB)(la) ~ 27 mag arcsec, equivalent to a stellar mass surface density of ~ 1 M pc -2 . 
S 4 G thus provides an unprecedented data set for the study of the distribution of mass and stellar 
structures in the local Universe. This large, unbiased and extremely deep sample of all Hubble types 
from dwarfs to spirals to ellipticals will allow for detailed structural studies, not only as a function 
of stellar mass, but also as a function of the local environment. The data from this survey will serve 
as a vital testbed for cosmological simulations predicting the stellar mass properties of present-day 
galaxies. This paper introduces the survey, describes the sample selection, the significance of the 
3.6 and 4.5/im bands for this study, and the data collection & survey strategy. We describe the 
S 4 G data analysis pipeline and present measurements for a first set of galaxies, observed in both the 
cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, 
position angle, axial ratio, inclination at //3.6 J um(^-B)= 25.5 and 26.5 mag arcsec -2 (equivalent to 
« /xb(^4-B)=27.2 and 28.2 mag arcsec -2 , respectively). These measurements will form the initial S 4 G 
catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged 
radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each 
galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, 
and nuclear point source, where necessary. Together these data products will provide a comprehensive 
and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe and 
enable a variety of scientific investigations, some of which are highlighted in this introductory S 4 G 
survey paper. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies:structure 
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1. INTRODUCTION 

Understanding the distribution of stars within a galaxy 
is akin to the study of the endoskeleton of a body; embed- 
ded within the galaxy is the fossil record of the assembly 
history and evolutionary processes of cosmic time. The 
first step in unraveling this history is to obtain a com- 
plete census of the stellar structures in galaxies in the 
local volume. This is the primary motivation for the 
Spitzer Survey of Stellar Structures in Galaxies (S 4 G). 
S 4 G is a volume- limited (d < 40 Mpc), size-limited 
(D25 >!') and apparent B-band brightness (corrected 
for inclination, galactic extinction and K-correction) lim- 
ited (niB rarr <15.5) survey of 2,331 nearby galaxies at 
3.6 and 4.5/xm wi th the Infrared Array Camera (IRAC) 
Fazio et aLll2004D o n the Spitzer Space Telescope (SSC) 



(jWerner et alJ feOCM). 



Over the last 70 years, numerous surveys have sought 
to create baseline data sets for n earby galaxies with as 
few as t ens of galaxies (ANGS T: iDalcanton et al.l l2009t 
SINGS, iKennicutt et all 12003ft to tens of thousands of 
galaxies (e.g.. RC3. Ide Vaucouleurs et al.l 119911) . While 
these surveys was designed specifically for particular sci- 
entific goals, none of them provides an accurate inven- 
tory of the stellar mass and the stellar structures in 
nearby galaxies. The main reason is that infrared sur- 
veys, where light from the old stellar population is bet- 
ter measured due to reduced dust extinction and con- 
tamination from star formation, are extremely difficult 
to conduct from the ground. The largest ground-b ased 
infra red surveys (e.g., 2MASS.lSkrutskie et al.ll2006b DE- 
NIS, llcchteinieLaUlilll iPaturel et alll2003l ) are rather 
shallow in depth, and as a result, studies of nearby galax- 
ies have been restricted to high surface brightness inner 
disks. In contrast, deeper infrared surveys (e.g., OS- 
UBGS. lEskridge et alll200l NIRSOS. lLaurikainen et all 
l20inh have imaged only a few hundred galaxies. S 4 G is 
specifically designed to provide over an order of magni- 
tude improvement in the sample size and several magni- 
tudes deeper data than existing surveys. 

S 4 G builds upon the two previ ous Spitzer Legacy sur- 
veys of nearby galaxies: SINGS (jKennicutt et all 20031) 
and t he Local Volume Legacy Survey (LVL) (|Lee et all 
2008). SINGS was designed to probe star formation, 
dust and polycyclic aromatic hydrocarbon (PAH) emis- 
sion in all representative environments within galaxies in 
a sample of 75 objects. The LVL survey was designed 
to study spatially-resolved star formation and the red 
stellar population within galaxies in a local volume of 
11 Mpc. The local volume limits the 258 LVL galaxies 
to consist mostly of dwarf, irregular and late type sys- 
tems (see green histogram in Fig. 1). The remainder of 
the Spitzer archive from the cryogenic mission is inade- 
quate for the science goals outlined here because it lacks 
sufficient numbers of galaxies of all stellar masses, partic- 
ularly at M* < 10 9 Mq, Hubble types and environments. 
As shown in Figures [TJ and [3J S 4 G explores the full range 
of stellar structures in a representative and large sample 
of galaxies of all types, masses and in diverse environ- 
ments. S 4 G is designed to provide the ultimate baseline 
data set for the study of stellar structures and mass in 
nearby galaxies. 

In this paper we introduce S G, describe the data col- 
lection and analysis strategy and briefly describe the 




Fig. 1. — Left: Histogram of all galaxies in S 4 G as a func- 
tion of the Hubble T-type shown in black. The green and blue 
filled histograms show archival data from the LVL (green) and 
other GO/GTO data (blue). Right: Histogram of masses for the 
S 4 G sample galaxies shown in black, where masses are calculated 
from the 2MAS S photometry assuming the M/L function from 
(IBell et alJ200a) . The biggest gains in building a statistically com- 
plete sample are for galaxies with log (M/Mq) < 9.2. 




Fig. 2.— Galaxies from S 4 G plotted on the lo cal large scale 
structure map from the 2MASS XSC ( Jarrctt 2004) shown in pur- 
ple. The red circles show the warm mission S 4 G galaxies to be 
observed. The cyan and yellow circles show LVL and GO/GTO 
archival galaxies. With the full S 4 G sample, we will have the abil- 
ity to make a less biased study of the galaxy host properties as a 
function of the local large scale structures. 

variety of scientific investigations possible with these 
data. All S 4 G galaxies are processed through a uniform 
pipeline to create the deepest and largest mid-infrared 
image catalog of nearby galaxies to date. For every 
galaxy we measure the magnitude and surface bright- 
ness and tabulate the galaxy diameter, position angle, 
axial ratio, and inclination at ns.^^m (AB) = 25.5 and 
26.5 mag arcsec~ 2 . We also compute the azimuthally- 
averaged radial profiles of ellipticity, position angle, sur- 
face brightness and color. This fundamental catalog of 
galaxy properties will allow us to classify their mor- 
phology based on their ste llar structu r es - a prelimi- 
nary study is presented in iButa et "all (|2010f) . In the 
long term, the main goal of the S 4 G project is to ob- 
tain a detailed understanding of the properties of the 
different stellar structures, their formation and evolu- 
tionary paths and their role in the broader picture of 
galaxy evolution. This will be done by decomposing 
every galaxy into its constituent structural components 
(e.g., bulge /sp heroid, bar, disk, nucle ar point source) us- 
ing GALFIT (IPeng et al.M200a l20loh . 

A number of scientific studies are already under way 
by the S 4 G team. These include a study of outer disks 
to determine the variety and frequency of different disk 
profile(s), truncations or thresholds and their relation- 
ship to the build up of galaxy disks, a study of shells 
and debris around galaxies and their relationship to 
central star formation or AGN activity, an analysis of 
the structural properties of bulges/spheroids (classical, 
disky, boxy /peanut-shaped), determination of the bar 
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fraction and bar properties in relationship to the host 
galaxies and environment, an analysis of the stellar arm- 
interarm variations in spirals, an investigation of the un- 
derlying old stellar population in rings and a search for 
fossil rings, an analysis of the old stellar content of dwarf 
galaxies and its implication for their formation history, 
etc. Longer term studies include star formation history 
and assembly of galaxies using GALEX and Ha data, 
the Tully- Fisher relationship usin g existing HI data (e.g., 
ALFALFA. rGiovanelli et al.ll2005L HIPASS. lBarnes et al.1 
20011 ). and testing of new mid-infrared diagnostics of 
AGN. We describe these studies in more detail in 
Ultimately, S 4 G will lead to numerous other studies of 
astrophysical phenomena in the nearby Universe and pro- 
vide a key anchoring data set for studies of galaxy evo- 
lution. 

2. THE SPITZER 3.6 AND 4.5/jM BANDS: BEST TRACERS 
OF STELLAR MASS 

The 3.6 and 4.5/xm IRAC bands on Spitzer are ideal 
tracers of the stellar mass distribution in galaxies because 
they image the Rayleigh- Jeans limit of the blackbody 
emission for stars with T > 2000 K. Moreover, the mid- 
infrared [3. 6] -[4. 5] color in nearby galaxies is nearly con- 
stant with radius, and independent o f the age of the ste l- 
lar population and its mass function (P ahre et al.l f2004V 
Analysis of the tilt in the fundamental plane for ellipti- 
cals shows that the Spitzer IRAC bands are better tracers 
of the stellar mass than even K-band emission, which is 
more sensitive to variations in the underlying stellar pop- 
ulation; the IRAC data are more uniform and therefore 
are th e best and perh aps most direct tracers of stellar 
mass (jJun fc Iml [2008) . Additionally, on Spitzer these 
bands are so sensitive that S 4 G will image to extremely 
low stellar mass surface densities (~ IM0PC -2 ). 

There are minor contaminants to the old stellar pop- 
ulation light in the 3.6 and 4.5^m IRAC bands. In the 
3.6/im band, there is a very weak 3.3/im PAH feature, 
which contributes negligibly to the overall emission in 
the band (< 2% as the equivalent width is ~0.02^m, 
iTokunaga et~aHll991[ L Hot dust (T d >500K) from very 
small grains can contribute light to both bands with a 
higher contribution in the 4.5/im band, but this condi- 
tion occurs only near active galactic nuclei or extreme 
starbursts. By examining the colors of the galaxies, we 
should be able to use the two IRAC bands to remove the 
effects of the 3.3/im PAH or hot dust emission. Another 
important advantage of having 2 bands is that by co- 
adding them, we can further improve the signal to noise 
and push the study to even fainter regions of the galax- 
ies. In both bands, the contamination from young red 
supergiants is also very low. From previous studies, it 
is known that in the near-infrared, the overall contam- 
ination from young red supergiants to a galaxy is only 
~3% but in regions of st ar formation, the contribution 
may be as high as 25% (|RhoadslH9M ). These regions, 
however, are easily identifiable and can be flagged. At 
the IRAC wavelengths, the fractional contribution from 
star-forming regions is expected to be the same (or lower) 
because these wavebands are at longer wavelengths on 
the Rayleigh- Jeans tail of the spectral energy distribu- 
tion. At these wavele ngths, there is lower dust extinction 
(|Draine fc Led [T984D and no other significant emission 
sources. Thus, with the extremely low surface bright- 



ness limits, S G offers a unique and virtually dust-free 
view of the distribution of mass in stellar structure in the 
nearby Universe. 

Although ground-based near-infrared observations can 
offer higher angular resolution than Spitzer, the main 
obstacle for ground-based observations is the very high 
and variable sky brightness (typically [Ik ~ 13.4 mag 
arcsec -2 ). The surface brightness level in the IRAC data 
is over 10 magnitudes below the typical sky brightness 
level. To get to this level from the ground, one would 
need to characterize the sky brightness, the flat field and 
instrumental variations to better than 0.0009%. This 
is currently not possible with any existing (or planned) 
near-infrared survey. Surveys like UKIDSS LAS or 
VISTA VHS are planning to reach a depth of 18.4 (21.2 
in AB) mag arcsec -2 in the K s band, which translates 
to 30 L0pc~ 2 (even with the brighter stellar emission at 
K s band), whereas S 4 G reaches ~ 2 L & pc -2 . It is the 
stability of the background that is the biggest advantage 
for space-based observations. 

3. SAMPLE SELECTION 

We chose all galaxies with radial velocity V radio < 3000 
km/s (corresponding to a distance d < 40 Mpc for a 
Hubble constant of 75 km/s/Mpc), total corrected blue 
magnitude niB corr < 15.5, blue light isophotal angular 
diameter D25 > 1.0 arcmin at galactic la titude |b| > 30° 
using HyperLEDA (jPaturel et ; all 12003( 1 . The choice of 
V ra rfi does limit the sample to galaxies with a radio- 
derived (e.g., HI) radial velocities in HyperLEDA, and 
as a result misses some galaxies for which only optically- 
derived radial velocities exists; currently, a sample cho- 
sen with V optical instead of V ra dio would contain 2997 
galaxies. A comparison of the galaxy properties of the 
two samples shows that the S 4 G sample misses galaxies 
that are small, relatively faint and early-type (gas-poor) 
systems from the volume surveyed. 

We also note that our sample was defined using Hy- 
perLEDA in September 2007 - since then 98 more galax- 
ies have been added to the HyperLEDA data base that 
meet the S 4 G criteria described above and we expect that 
some more may be added as better data become avail- 
able. Given the limited lifetime of Spitzer, we are unable 
to go back and acquire data on these additional galaxies 
but we expect that the few percent additional galaxies 
will not strongly influence the core characteristics of the 
sample, except for the bias against early-type galaxies as 
noted above. 

The choice of a 40 Mpc volume is arbitrary - it was 
chosen to be large enough to provide a statistically sig- 
nificant number of galaxies of all types and to be repre- 
sentative of a large range of the local large scale structure 
environment (as shown in Fig. [2]). Our exp erience with 
previ ous surveys of galaxy properties (e.g., iSheth et al.l 
2008) has shown that a few thousand galaxies are needed 
to accurately measure and account for completeness ef- 
fects in mass, color and size selection of galaxy samples 
for robust investigations of relationships between galac- 
tic structure and host galaxy properties. The size cut, 
log D25 > 1.0, was made to ensure that galaxies were 
large enough for a detailed study of their internal struc- 
ture (at 40 Mpc, 1' ^11.6 kpc). The apparent size and 
apparent magnitude cut were chosen to match the RC3 
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limits. The galactic latitude cut, |b| > 30°, was used 
to minimize the unresolved Galactic light contribution 
from the Milky Way disk. The full sample of 2,331 S 4 G 
galaxies is listed in Table [1] In Figure Q] we show the 
distribution of S 4 G galaxies as a function of Hubble type 
and stellar mass, along with the existing distribution of 
the LVL galaxies (green) and GO/GTO galaxies (blue) 
including the SINGS galaxies. There are some notewor- 
thy sub-samples within S 4 G . There are 188 early-type 
galaxies (ellipticals and lenticulars), 206 dwarf galaxies 
(defined as Mb > -17), of which 135 are DDO dwarfs and 
465 edge-on (i > 75°) systems. Also over 800 of the S 4 G 
galaxies are mapped beyond a diameter of 3XD25. These 
sub-samples provide the deepest, largest and most homo- 
geneous samples for a multitude of galaxy evolutionary 
studies, some of which are described in [JB] 

3.1. Archival / Cryogenic Spitzer Data 

Of the 2,331 galaxies in the S 4 G sample, 597 had some 
existing data at 3.6 and 4.5/mr from the Spitzer cryogenic 
mission. 125 of these are part of the LVL survey and 56 
are from the SINGS survey, both of which used the same 
observing strategy (Q we employ for S 4 G . Almost all of 
the archival maps have at least 240 s of integration time 
per pixel and are sufficiently deep. The only exceptions 
are NGC 5457 (96s), NGC 0470, NGC 0474 (150s), and 
NGC 5218, NGC 5216, and NGC 5576 (192s). 

In terms of area, 82 of the archival galaxies were 
mapped between 1.0 and I.5XD25, and 43 others were 
mapped to <1.0 XD25. While outer disk science for these 
125 galaxies would benefit from an extended map, they 
represent only a small fraction of the total sample, and 
the repetition of observations would have yielded only 
an incremental scientific return. Therefore, we decided 
not to repeat observations of any of the archival galax- 
ies. We also decided not to map the SMC and the LMC, 
which meet our selection criteria but are very large and 
are being observed as part of other GO programs. M33, 
the other commonly studied Local Group galaxy was ob- 
served during the cryogenic mission and is included in 
the archival part of our sample. All other galaxies that 
were found via the HyperLEDA data base search in 2007 
September are observed (see note earlier about the mod- 
est additions to the HyperLEDA database in the last 
three years). Although we do not repeat observations of 
any of the archival galaxies, we are processing all galaxies 
in the S 4 G sample, except the SMC and LMC, through 
the same S 4 G reduction and analysis pipelines described 
in gi 

4. SPITZER POST-CRYOGENIC MISSION OBSERVATIONS 
STRATEGY 

The 1,734 S 4 G galaxies being observed in the post- 
cryogenic mission are mapped using either a small 
dithered map or a mosaicked observation with both 
IRAC channels. All galaxies are observed with a total 
on-source integration time of 240 s which leads to an rms 
noise of /i(lo-) - 0.0072 and 0.0093 MJy sr" 1 at 3.6 and 
4.5/im, respectively. This translates to a typical surface 
brightness sensitivity of hab ~27 mag arcsec -2 or a stel- 
lar surfa ce density of 1. 13 M©pc -2 (assuming a solar 
= 3.33. IWOTthev([T99l M/L^=l. iHeraudeau fc Shmenl 
H997I) . 




XrM< 



Fig. 3.— Left: A typical dithering AOR shown for NGC 2919 
in the warm mission. The blue (red) overlay shows the 3.6/^m 
(4.5 fim) coverage. The galaxy is mapped in two separate AORs 
separated by at least 30 days. The flanking fields are reduced by 
the S 4 G pipeline and may also contain serendipitous objects such 
as interacting galaxies or tidal debris. Right: A typical mosaicking 
AOR shown for NGC 1515 in the warm mission. The blue grid 
shows the coverage of the 3.6/im camera and the red grid shows the 
coverage of the 4.5 jira camera. The galaxy is mapped to >1.5xD25 
by both arrays. The mosaicking pattern takes steps of 146.6" with 
30s integrations in each location. The total 240 s integration time 
is obtained in two separate AORs. Note that significant amounts 
of sky are also mapped on opposite ends of the galaxy by each of 
the individual arrays. 

1560 of the new observations are of galaxies with D25 < 
3.3'. We can map these galaxies to I.5XD25 with a 
single pointing because the Spitzer field of view is 5'. 
Hence all of these galaxies are mapped using a stan- 
dard cycling small dither pattern with 4 exposures of 30s 
each in two separate Astronomical Observation Requests 
(AORs). Each pair of AORs is separated by at least 30 
days to allow for sufficient rotation of the telescope so 
that the galaxy is imaged at two distinct orientations. 
Two AORs allow us to better remove cosmic rays, and 
the redundant information gathered by the two visits al- 
lows us to characterize and remove image artifacts (e.g., 
muxbleed, column pulldown, etc (see the IRAC Instru- 
ment HandboolP^f) and possible asteroids. Drizzling the 
data from the two visits also allows us to achieve sub- 
pixel sampling with which we can reconstruct images 
with better fidelity than would be possible from a sin- 
gle visit. We also note that since both the 3.6 and 4.5^m 
arrays collect data simultaneously, their offset placement 
in the Spitzer focal plane means that an additional flank- 
ing map is made adjacent to the galaxy in each of the two 
channels (in some cases, the flanking field from each of 
the AORs may not overlap exactly due to the telescope 
rotation, but this does not affect the observations of the 
main galaxy) . An example of a typical dithering AOR is 
shown in the left panel in Figure [3] 

One of the considerations in designing these AORs was 
the effect of saturation from bright stars on the chip and 
the scattered light from bright stars falling in one of the 
three scattering zones surrounding the arrays0 Stars 
brighter than 5th magnitude can leave a latent charge 
on the chip. Stars brighter than 11th magnitude in the 
scattering zones may scatter light into the main field of 
view of the IRAC chips. To mitigate these cases, we used 

23 See http : //ssc . spitzer . caltech . edu/irac/iracinstrumenthandbook/ 

24 A description of this problem is 
given in the IRAC instrument handbook at: 

http: / /ssc. spitzer. caltech.edu/irac/iracinstrumenthandbook/ 
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the Spitzer Observing Tool (SPOTS) to examine every 
galaxy with a star, uyks < 8 (stars fainter than vuks 
— 8 do not contribute significantly to the background) 
within an area of 300", covering the IRAC field of view 
and scattering zones adequately, for all possible visibility 
windows through 2011. In a majority of these cases, no 
modifications of AORs were needed because the offend- 
ing star(s) fell out of the chip or the scattering zones. 
In the remaining cases, we constructed AORs with tim- 
ing constraints, slightly offset pointings or used medium 
dithers or mosaics to reduce the effect of the bright stars. 

For galaxies with D25 > 3.3', we mosaic them in array 
coordinates with offsets of 146.6" with 30s integrations 
at each location to create a map >1.5xD25. This leads 
to a map with each pixel observed four times at all wave- 
lengths in the core of the map. In addition, on each side 
of the galaxy there will be a region where each pixel is ob- 
served twice but only in one of the two channels. Like the 
dithered AORs, each mosaic is observed twice but with 
only a follow-on constraint so that each mosaic overlaps 
the other closely. An example of this type of an AOR is 
shown in the right panel of Figure [3] 

Although larger mosaics could be made, our expe- 
rience with SINGS data showed that the dominant 
noise term comes from variations of the background 
(jRegan fc Teubenl [200l iRegan et all [2006h . The best 
background for a galaxy includes sufficient sky not far 
from the area of interest. Our map sizes are thus fine- 
tuned to achieve the best possible signal-to-noise for all 
galaxies. 

5. THE S 4 G PIPELINE 

The S 4 G pipeline is divided into four distinct parts, 
which we refer to as Pipelines 1, 2, 3 and 4. Pipeline 1 
takes the raw basic calibrated data (BCD) fits files and 
converts them to science-ready data. Pipeline 2 creates a 
mask for the foreground and background objects for each 
of the individual galaxies. Pipeline 3 measures the stan- 
dard global galaxy properties such as size, axial ratio, 
magnitude, color, etc. and computes the radial profiles 
of the standard properties. Pipeline 4 deconstructs each 
galaxy into its major structural components. Each of 
these is described in detail here. 

5.1. S^G Pipeline 1 

The purpose of pipeline 1 (PI) is to create science- 
ready mosaics from the two visits of each target. To 
do this PI performs two major steps. First, it matches 
the background level in the individual images to account 
for drifts in the zero point of the amplifiers. Second, 
it creates the final mosaic using the Drizzle package of 
IRAI0 

To match the background levels in the individual im- 
ages we find the regions of overlap between all pairs 
of overlapping images. Typically we require that the 
overlap region contain 20000 pixels. For all of the new 
post-cryogenic mission observations (and archival obser- 
vations that used the same observing strategies), such 

25 http: //ssc . spitzer . caltech. edu/warmmission/propkit/spot/ 

26 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



an overlap is always available. However, when reducing 
other archival data with less overlap between the frames, 
we reduce the number of pixels in the overlap region to 
as low as 5000 pixels - this still provides an adequate 
number of pixels for estimating the background. 

Within these regions we determine the brightness level 
of the 20th percentile pixel. We use the 20th percentile 
brightness level to better avoid detector artifacts such 
as muxbleed which can affect even the median. By us- 
ing the 20th percentile brightness level we have a better 
chance of using the true sky level. We then perform 
a least square solution using the brightness differences 
between all the pairs of overlapping images, minimizing 
the residual background difference. Since we use the dif- 
ference, we need to use the background level in the first 
image to set the zero point of the solution. We then make 
corrected images by adding the solved-for corrections and 
use these images in the formation of the final mosaic. We 
create the final mosaic following the standard prescrip- 
tion of the STSDAS DITHER package. This method 
removes the cosmic rays by first forming a mosaic from a 
median combination of the images and then comparing 
the individual images to the median mosaic and flagging 
any cosmic rays. The final mosaic is formed by drizzling 
the individual images. The resulting mosaic has a pixel 
scale of 0.75" and we correct for the change in pixel size 
to keep the units of MJy sr _1 . 

The relative astrometry of the mosaicked images is ex- 
cellent because the Spitzer pipeline already incorporates 
the 2MASS positions to update its astrometry for each 
of the basic calibrated data frames. The S 4 G pipeline 
therefore does not require the type of cross-correlation 
that was needed previously for the SINGS pipeline to 
align overlapping fields. The relative astrometric accu- 
racy is <0'.'l. 

The point spread function for IRAC is asymmetric and 
depends on the spacecraft orientation angle. Pipeline 1 
images are produced by combining two different visits to 
the galaxy and therefore we estimate the PSF by com- 
bining stars in the foreground and background for six 
typical S 4 G galaxies to create a "super"-PSF which has 
a typical FWHM of 1"7 and 1"6 at 3.6/jm and 4.5/xm 
respectively. For the most distant galaxies in our survey 
at a distance of ~40Mpc, this resolution corresponds to 
a linear scale of ~300pc. At the median distance of the 
galaxies in this survey at 21.6 Mpc, this corresponds to 
a linear scale of ~170pc. 

5.2. S 4 G Pipeline 2 

After the science-ready images are produced by the 
S 4 G pipelin e, we generate masks for point sources using 
SExtractor (jBertin fe Arnoutsl [T996) for both channels. 
Each mask is checked by eye and iterated to identify any 
point source otherwise missed by SExtractor. We also 
unmask any region of the galaxy that may be have been 
incorrectly identified by SExtractor. The images with 
their corresponding masks are then run through the third 
part of the S 4 G pipeline (Pipeline 3, P3 hereafter). 

5.3. S^G Pipeline 3 

The first step in P3 is to determine the local sky level 
around each galaxy. We want to be sufficiently far from 
the galaxy to ensure that there is no contamination from 
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the galaxy but not so far that the background is not 
truly local. We do this by computing the median sky 
value in two concentric elliptical annuli centered at the 
position of the galaxy. Typically we start the inner an- 
nulus at a distance of I.5XD25. This annulus is divided 
in azimuth into 45 regions (or sky boxes). Each of these 
regions is then grown outwards until it encompasses a 
total of 4000 non- masked pixels^ 7 !, as shown in Figure |4j 
The second annulus begins at the end of the first annu- 
lus and the same process is repeated. We compare the 
measured sky values in the two annuli for any radial gra- 
dients which might be indicative of light contamination 
from the galaxy. If we find that there is contamination, 
then we move the inner radius outwards and repeat the 
process until we are assured that we are sampling the 
local sky. 

In our analysis so far, the inner radius of the inner 
annulus has varied between 1.2 and 2XD25. In all cases 
for the warm mission data and a majority of the archival 
data, it has been sufficent to set the inner radius of the 
annulus to either 1.5 x or 2x D25. However in a few 
cases where the archival data are shallow and the galaxy 
is relatively blue, we moved the inner radius to 1.2 XD25 
to be immediately outside the galaxy to get as accurate a 
measurement possible for the local sky level. Automating 
this procedure is non-trivial because variations at these 
faint levels (/^3.6/jm >27 (AB) mag arcsec -2 ) can be from 
a lopsided galaxy, debris or tidal structures, or variation 
in the background. Therefore we examine every galaxy 
by eye to ensure that contamination from the galaxy light 
to the measurement of the sky background is as low as 
possible and move the radii outwards as needed. 

The sky background computation method for three of 
the galaxies in the S 4 G sample (NGC 0337, NGC 4450, 
NGC 4579) is shown in Figure |U In the central column 
of this figure, we show an example case for NGC 4450 
where the median sky in some of the boxes is affected 
by the bleeding of flux from a (masked) but saturated 
nearby bright field star. In such cases we ignore the af- 
fected pixels or move the annuli further out to get the 
best estimate for the background. Finally, we use the 
standard deviation obtained from the 90 sky boxes with 
the average standard deviation within each box to esti- 
mate the contribution of the low-frequency flat-fielding 
errors and possible additions to the total error budget in 
the background determination. 

The typical sky brightness level is fainter than 26 
(AB) mag arcsec -2 at 3.6 and 4.5 fim. The Poisson noise 
from the galaxy dominates the total error budget in the 
central regions, while errors in the sky background are 
the main source of uncertainty in the outer parts. Large- 
scale background errors can be particularly significant 
in galaxies with large apparent sizes. However, in the 
fainter outer parts, the noise is significantly reduced from 
azimuthally averaging over a large number of pixels in 
the outer galaxy - we are therefore able to measure the 
galaxy profiles well below the typical sky brightness. The 
typical S/N for a S 4 G galaxy at the canonical ^3.6^ = 
26.5 (AB) mag arcsec -2 is > 3. 

The images resulting from Pipeline 1 are calibrated in 
units of MJy sr -1 . The conversion from MJy sr -1 to AB 

27 For archival galaxies that fill the mosaicked frame, this re- 
quirement is reduced to 1000 non-masked pixels. 
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Fig. 4. — Results of the determination of the sky background 
for the S 4 G galaxy NGC 0337, NGC 4450 and NGC4579. Top 
row: Results for the IRAC 3.6/mi image for each galaxy. The D25 
ellipse is shown in red and inner edge of the annulus from where 
we begin to estimate the sky is shown in green. Second row: The 
elliptical annuli used to compute the sky background are shown. 
The red segments indicate the 45 regions of the inner annulus. The 
outer annulus used to trace the gradient is shown shaded by the 
grey scale annulus. The changing shade from light grey to black 
corresponds to the 45 regions of the outer annulus. In this row all 
the masked pixels are shown in black but note that the black color 
of the outer annulus does not mean that these regions arc masked. 
Third row: The median and standard deviation for each of the 90 
sky boxes is shown here. Bottom row: A histograms of the pixel 
values within each of the 90 sky boxes is shown in this panel. 

magnitudes is such that the zero-point to convert fl uxes 
in Jy to magnitudes does not vary with wavelength (|Okd 
[1971 : 

77i AB (mag) = -2.51ogF„(Jy) + 8.9 (1) 

From this definition one can derive the corresponding 
expression for surface brightness: 

^AB(mag) = -2.5 log I v (MJy str -1 ) + 20.472 (2) 

5.3.1. PS: Measuring the Galaxy Host Properties 

After the sky level has been properly measured, we use 
the IRAF routine ellipse to determine the radial pro- 
files of intensity, surface brightness, ellipticity and po- 
sition angle for all galaxies in both bands. The center 
is determined using imcenter and kept fixed during the 
fit, whereas the ellipticity and position angle are left as 
free parameters. We generate profiles with two differ- 
ent radial resolutions, by incrementing the semi-major 
axis in 2" and 6" at each step, respectively. The lat- 
ter profiles have a coarser resolution and a better S/N 
ratio - these are used to measure the RC3-like parame- 
ters. The low- and high-frequency sky background errors 
described above, together with the Poisson noise of the 
source along each isophote, are considered together in 
the final er ror budget of the surface photometry at each 
radiu s (see iGil de Paz et al.l 120051 : iMunoz-Mateos et al.l 
2009). The errors in the position angle and ellipticity at 
each radius are computed by the ellipse task from the 
mis of the pixel values along each fitted isophote. The 
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Fig. 5. — Surface photometry results for the S 4 G galaxy 
NGC 0337. Top: Surface brightness profiles for 3.6^m (black filled 
circles) and 4.5/im+l mag (grey filled circles). The vertical dashed 
line shows the position of the B-band D25 major axis radius. Mid- 
dle: Variation of the ellipticity of the best-fitting isophote with 
the radius. For the sake of clarity we have applied an offset of 0.1 
to the 4.5/xm profile. Bottom: Position Angle of the best-fitting 
isophote as a function of radius. Note the 20 dg offset applied to 
the 4.5/im profile. The boxes drawn at the position of the D25 
major axis radius in the latter two panels indicate the ellipticity 
and PA listed in the RC3 catalog, respectively. 



corresponding errors for the values of these geometrical 
parameters at 25.5 and 26.5 (AB) magarcsec -2 are inter- 
polated between the corresponding adjacent values. The 
error in the semi-major axis is computed from the change 
in the radius affected by moving the surface brightness 
profile by the measured errors and an additional 10% 
error added in quadrature to the measured errors. The 
10% error is an estimate of the aperture correction uncer- 
tainty as estimated by the Spitzer Science Center instru- 
ment team. We compute the error in the semi-major axis 
for each galaxy. Sample results for NGC 0337 and NGC 
4579 are noted below. For NGC 4450, the scattered light 
from a background star strongly affects the computation 
of the surface brightness error and therefore the error in 
the semi-major axis that is unrealistic. In such cases, the 
error in the semi-major axis can be assumed to be the 
median semi-major axis error in the S G sample. The 
maximum error in the determination of the semi-major 
axis is 6", the step size of the coarse resolution fitting pro- 
cedure. In Figures El E] and [7] we show three sample sets 
of profiles for three galaxies in the sample: NGC 0337 
(SBd), NGC 4450 (SAab) and NGC 4579 (SABb). The 
basic RC3-like galaxy properties (semi-major axis, axial 
ratio and position angle) are tabulated at the levels of 
25.5 and 26.5 (AB) mag arcsec -2 in both the 3.6^m 
and the 4.5/im images, and are quoted in Table [2] 
Along with the surface photometry we also measure 
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The same as Figured] for NGC 4450. 
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Fig. 7. — The same as Figure [5] for NGC 4579. 
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the curve of growth (e.g., iMunoz- Mateos et all :2009; 
iGil de Paz et al.l 120071 ) with the integrated flux up to 
each radius. By fitting the accumulated magnitude as 
a function of the magnitude gradient at each point we 
obtain the asymptotic magnitude as the y-intercept of 
that fit (see IGil de Paz et~al]|2007|) . Once the total mag- 
nitude is known, it is straightforward to locate the radii 
containing a given percentage of the total galaxy hum- 
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Fig. 8. — Five-component decomposition model for NGC 1097. 
The panels show the observed 3.6 micron image (left), the 
model image (middle), and the normalized residual image (OBS- 
MODEL) /MODEL (right). North is up and east is to the left per 
the usual conventions and the units are in arcseconds. The decom- 
position includes bulge, disk, bar, central point source and nuclear 
ring components, as illustrated in Fig [9] Note how the residual 
image emphasizes the outer spiral and the inner ring structures. 

NGC 1097 

121 

BULGE 0.200 

14 TOT MODEL 

BAR 0.212 




1 10 100 1000 

Distance along sky (arcsec) 



Fig. 9. — The different sub-components of the NGC1097 fit in 
Fig. [8] The white dots show the magnitude of each pixel of the 
observed image, as a function of distance from the galaxy cen- 
ter along the sky plane, while the dark points indicate the total 
model. The different colors indicate different sub-components: the 
numbers after the labels indicate the relative contribution of the 
component to the total model flux. In this type of plot axially 
symmetric components appear as curves, while non-axisymmetric 
components (i.e. the bar) fill a region confined by their major and 
minor axis profiles. Similar plots ha ve been used t o visualize BD - 
BAR l ILaurikainen et al.ll2005l , I2009T ) and BUDDA HGadottill20"08h 
decompositions. 

nosity, from which concentration indices can be deter- 
mined. In Table [3] we quote the asymptotic magnitudes 
of the three sample galaxies considered here, to gether 
with the C 3 i (|de Vaucouleurs et alii 19 771 ) and Cap (|KentJ 
1985) concentration indices. 

5.4. S^G Pipeline 4 

In Pipeline 4, we decompose the two-dimensional 
stellar distribution in each galaxy into di fferent sub- 
component s using the version 3 of GALFIT (|Peng et al.l 
I2002L 120101 ) . GALFIT is a parametric fitting algorithm 
that offers a large flexibility on the fitted models; the op- 
timal solution for the model parameters is found with the 
Levenberg-Marquardt algorithm, performing non-linear 
least-squares % 2 minimization of the difference between 
observed and model image s. For our pip eline, we could 
have also used BUDDA . (|Gadottil 120081 ) and BDBAR 
(jLaurikainen et aT1l2005l ) . We tested all three algorithms 
using a set of galaxy images and verified that GALFIT 
produces the same results as BUDDA and BDBAR if the 
three codes are run with the same set of input images, 
models and initial parameters. We chose GALFIT be- 
cause it is currently the most sophisticated algorithm, 
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NGC 1 0&7 model with spirals 




Fig. 10.— Illustration of the capabilities of GALFIT3.0: a de- 
composition for NGC 1097 including spiral sub-components. While 
for many scientific purposes it is reasonable to use parametric func- 
tions like Sersic or Ferrers functions for fitting structural compo- 
nents, GALFIT also allows modeling of curved or irregular struc- 
tures. As an illustration of modeling where the spiral arms form 
part of the decomposition, a model for NGC 1097 is shown - it 
shows a remarkable resemblance to the observed image. We cau- 
tion the reader however that since the residuals contain contribu- 
tions from partially modeled spiral arms, the model parameters 
output by GALFIT are not necessarily ideal for the study of faint 
structures in galaxies. 

and also well known by a large portion of the astronom- 
ical community Other studies in whi ch softw a re wa s 
developed for s i milar purposes include iSimardl (119981) : 
iPignatelli etaLl (|2006D : iMendez-Abreu et al.l (|2008D ~ 

To simplify the use of GALFIT for a large num- 
ber of cases, a set of new IDL-based tools have been 
created ("GALFIDL"). They include automatic cre- 
ation of the input files for decompositions, contain- 
ing reasonable first guesses for the initial parameters 
based on P3 products, as well as routines for run- 
ning and visualizing the decompositions. GALFIDL 
will be available to the community concurrently with 
the publication of the data and the GALFIT analyses; 
an overview of these p rocedures is already available at 
|http : //cc . oulu . f i/ ^ hsalo/galf idl . html. 

As a standard part of the pipeline analysis, decomposi- 
tions include bulges and disks, and, if appropriate, bars, 
nuclear components, and in some cases multiple exponen- 
tial disks. Taking into account that decompositions, par- 
ticularly at high redshift, are often made in a more simple 
manner, we also compute the 2-component bulge/disk 
decompositions, and 1-component Sersic function fits. A 
decomposition for NGC 1097 is shown in Figs. [8] and 
[HI the 5-component model includes an exponential disk, 
a Sersic function for the bulge, and a Ferrers function 
for the bar. The galaxy has an extremely prominent nu- 
clear ring, which is fit using a Gaussian function, and a 
nuclear point source using the corresponding PSF. With- 
out deep, almost dust-free observations, like those in S 4 G 
, the measurement of bulge structural parameters, like 
shape, size and profile, might in some cases be compro- 
mised if the quantification of underlying structures, such 
as the disks, bars and rings is inaccurate. GALFIT3.0 al- 
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lows for even more sophisticated decompositions, which 
will be made for selected sub-samples of S 4 G galaxies. 
For example, it is possible to add more components, or 
one may want to fit the spiral arms, of which an example 
is shown in Fig. [TO] 

6. SCIENCE INVESTIGATIONS 

In this section we highlight some preliminary results 
from a number of science investigations that are and will 
be carried out by the S 4 G team. We also discuss possible 
studies that will be enabled by the S 4 G dataset. 

6.1. The Faint Outskirts of Galaxies 

A major area of study for S 4 G is a quantitative analysis 
of stellar structures in the faint outer regions of galaxies 
where the stellar surface density we trace is extremely 
low and comparable or lower than the atomic gas reser- 
voir. Already, deep optical studies have shown that, in a 
majority of spiral galaxies, the outer disks exhibit a sec- 
ondary, exponential compone nt. Stellar disks appear t o 
extend beyo n d 1.5 x D 2 5 (e.g.. iPohlen fc Truiillol (|2006f) : 
lErwin et al.l (|2008D ) with truncated o r anti-truncated 
profiles. GALEX (jMartin et al.1 120051 ) has found ex- 
tended UV (XUV) disks even farth er out, where it 



was ass umed only Hi gas ex isted ( e.g., [G il dc Pa z et al 
(|2005t) ; lThTlker et al.l (|2005l. I2007t) ; [Mufioz-Mateos et al 
(l200l : IZaritskv fc Ghristleinl (|200l . Not only is S 4 G 



significantly deeper than the previous optical data, but 
it also traces the stellar profiles into the H I disk of the 
galaxy - in over 800 S 4 G galaxies, we image an area 
>3D 25 . 

Another observation that has been much debated is 
the break between the inner and outer exponenti als pro- 
files (e.g.. lBakos eTani200l lAzzollini et al.ll2008l ). With 
S 4 G we are measuring the break radius and comparing 
it to the disk size and mass. The S 4 G catalog of these 
breaks will provide strong quantitative constraints for 
galaxy evolutionary models. In particular, this catalog 
will be key to quantify the possible role of radial stellar 
migration in shaping the outskirts o f disks ()Roskar et al.l 
120081 iSanchez-Blazquez et al.ll200l . 

In elliptical galaxies, deep optical observations have re- 
vealed the presence of fain t tidal tail and shell like struc- 
tures (e.g.. lBennert et al.l (|2008D ; iCanalizo et al.l (|2007l ). 
S 4 G contains 188 ellipticals and lenticulars - with this set 
and the S G depth we are quantifying the frequency and 
mass of shells and other debris features in the outskirts 
of ellipticals and lenticulars. 

6.2. Bulges: Classical, Disk-like, Boxy/Peanut 

Bulges are an inhomogeneous class of objects and 
differ ent types have been proposed in the litera - 
ture (IKormendv fc Kennicutd l 2004llAthanassoulall2005l ). 
IKormendv fc KennicuttI (12004 ) distinguish es bulges from 
pseudo-bulges, while lAthanassoulal (|2005l) distinguishes 
three categories: classical bulges, boxy /peanut bulges, 
which are a part of a bar, and disk-like bulges, 
formed out of disk material. Several studies have 
analyzed their observational prop erties and compared 
them with N-body model s (e.g. iKuiiken fc Merrifieldl 
19951: iCarqllo et all U997I: iBureau fc Freeman! I199S 



Gadotti 2009). Even in lenticular galaxies there seems 
to be evidence for disk-li ke bulges, based on struc- 
tural decompo s itions (e.g..lLaurikainen et alll2005l 120071 : 
iGadottil [200l iGraham fcWorlevI 120081) and kinemat - 
ics (e.g.. lPeletier et al.ll2007t IFalcon-Barroso et al.ll200ll . 
But until now few studies have analyzed the structural 
properties of bulges in a large number of galaxies; none 
have done it in a sample as large and deep as S 4 G . More- 
over, S 4 G encompasses a range of large scale structures, 
allowing us to examine the environmental influences on 
the formation of bulges. This may have important impli- 
cations since massive bulges are expected to form from 
mergers. 

From a prelim inary analysis of the S 4 G data, 
iButa et all (f2010h are finding that many late type ellip- 
ticals in the RC3 have disks, and thus should instead be 
classifie d as SOs, consistent with other deep near-infrared 
studies (|Laurikainen et aLll201Q[ ). With the S 4 G data we 
will examine how common the disk structure is in all 
types of ellipticals. 

Finally, the extremely deep data from S 4 G allows us 
to best determine the structural properties of the bulge 
component - in shallower data, estimates of shape, size 
and profile are uncertain when measurements of underly- 
ing components such as disks, bars and rings are impre- 
cise. As described in £15.41 the S 4 G project is carrying out 
detailed decomposition of all the galaxies, and expects to 
compile the most detailed data base of bulges of all kinds. 
The combined advantages of the S 4 G data (deep and in- 
sensitive to dust attenuation) and the uniform reduction 
pipeline (careful and detailed structural analysis) will 
allow for a suitable examination of galaxies with com- 
posite bulges, i.e., galaxies host i ng more than one bulge 
category (e.g.. | Erwin et al. 2003; Bar entine fc Kormendvl 
2009; No wak et all 120 101) . The spatial resolution of our 
survey is adequate for resolving typical bulges (our res- 
olution is ^100 pc at the median distance of the survey 
volume) but is inadequate for precise measurements of 
the bulge structure from the inability to resolve compact 
nuclear structures which are often present at the cen- 
ters of galaxies. These unresolved point sources affect 
the values of Sersic parameters, which on average are 
overestimated in grou nd-based optical images compared 
toHST studies (e.g., iCarollo et aT]|1997t iBalcells et all 
2003). We also expect that the reduced dust extinction 
may reveal triaxial structures in bulges, and thus allow 
us to probe beyond the standard symmetric Sersic fits. 

6.3. Bars 

Stellar bars are dynamically important struc- 
tures in galaxies. A bar-induced gas flow leads to 
a number of evolutionary effects from the trigger- 
ing of circumnuclear sta r bursts to the build - up o f 
bulges (IScoville fc HeTshl 119791 iSimkin et all I198L 
' Zaritskv fc Lolll986l: IAthanassoulalll992t IFriedli fe Benz 
1993; Martin fc Ro 
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Erwin et al.l 120031: IBureau fc Athanassoulall2005t iFisher 
20061: iDrorv fc Fished 120071: iFisher fc Drorvl 120081: 



The cosmological redshift evolution of the bar fraction 
is also an important signpost of the growth and dy- 
namic matu rity o f galaxies (e.g.. lAbraham et al.l 119991: 
Sheth et all 120031: fjogee et all [20041: lElmegreen et all 
20041: ISheth et all 120081) . Whereas the previous smaller 



and shallower studies gave mixed re sults for the c hange 
in the bar fraction with redshift, S heth et al.l (2008) 
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show that not only does the overall bar fraction in 
disk galaxies decline significantly over the last seven 
billion years, but also that the decline depends on 
the galaxy type. They show that the formation and 
evolution of bars is strongly correlated with the host 
galaxy mass, bulge-dominance and presumably their 
dark matter halo. The S 4 G data is critical to establish 
the precise local frequency of bars and its variation with 
galaxy type. It also provides a measu rement of the 
light fraction in the bar component (e.g.. IDurbala et al.l 
120081 : iGadottil 120081 : iWeinzirl et"aH 120091 ) and accurate 
measures of propertie s such as the bar length, shape 
and el liptici t y (e.g., iMenendez-Delmestre et all 120071 : 
lErwinl 120051 iGadottil l2010f) . Several techniques may 
be used to der ive bar strengths, b a sed o n calculating 
bar torques (jLaurikainen fc Salol l2002t IButa et all 
I2003T). and estimating the b ar-spiral arm contrast 
(jElmegreen fc ElmegreenI [1985). Analysis of the rel- 
ative Fourier intensity profiles of bars is compared 
with models of the evolution of barred galaxies (e.g. , 
lAthanassoula fe Misiriotisl l2002t iDebattista et all 120041) 
to constrain the history of th eir evolution and inte raction 
with their dark matter halo (j Athanassoulal [20031 ) . 

Another study of interest is the presence of sec- 
ondary or nuclear bars which may be responsible for 
feeding the central black holes and/or stellar clusters 
(Shi osman et~aT1[T990l ). The S 4 G observations, unaf- 
fected by the high and patchy extinction often present 
in galactic centers, allow us to detect all nuclear bars 
with sizes of r > 500 pc (our resolution, 2" ~ 200 pc at 
the median distance of S 4 G ) in galaxies of all types. 

6.4. Galactic Rings 

Outer rings are large, optically low surface brightness 
features that dominate the outer disks of some early- 
type barred and weak ly-barred galaxies (see review by 
IButa fe Combed ([19961 The properties of the rings con- 
strain their formation epoch, the dynamical time-scale 
at their respective radii, and the evolu tion of the bar 
pattern speed (e.g.. lAthanassoulal (|2003D . Finding outer 
rings around unbarred galaxies is also of great interest, 
because it would suggest that some unbarred galaxies 
may in fact be highly evolved former barred galaxies. 
Whereas ground-based near-IR imaging is rarely deep 
enough to detect outer rings reliably to study their stel- 
lar populations, colors, and other characteristics, the S 4 G 
data will not only allow us to detect many new outer 
rings, but also to make a complete census of these fea- 
tures in normal galaxies and to measure their stellar 
mass, shape, width and ellipticity. Likewise, S 4 G pro- 
vides an unprecedented view of circumnuclear and inner 
rings, which are typically found in barred spirals. These 
rings should have a broad, underlying older stellar pop- 
ulation background as they evolve, and it is this type of 
feature we will see in the S 4 G data. 

6.5. Spiral Arms 

A systematic study of spiral arm amplitudes as a func- 
tion of radius at 3.6 and 4.5/zm, free of confusing ef- 
fects from dust extinction, promises to initiate a new 
set of spiral arm studies. Although the discrepancies 
in the classical density wave m odel were solv ed with 
the swing amplification theory (jToomrel I1981D . there 
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Fig. 11. — Preliminary analysis of the arm-interarm ratios com- 
puted for a fiocculent spiral, NGC 7793, and a grand design spiral, 
M51, from the S 4 G 3.6/^m data. Note the increasing contrast in the 
grand design galaxy compared to the fiocculent galaxy. Detailed 
analysis of the stellar density variations will shed important light 
on spiral density wave theories. 
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Fig. 12. — Comparison of the M51 3.6 micron image with the 
N-body simulation model of Salo & Laurikaincn (2000a b). The 
right hand column shows the simulation at 40, 240 and 480 Myr 
from the top to the bottom panels respectively. In the model the 
near-IR, and rotation curve data existing at the time were used for 
setting up the initial disk and halo models. The simulation creates 
an outer southern tail when the companion (with mass ratio 0.55) 
crosses the primary disk towards the observer about 480 Myr ago. 
To obtain the observed apparent counter-rotating HI kinematics, 
the tail has a tilt of nearly 40 degrees with respect to the inner 
disk. The most recent disk plane crossing about 100 Myr ago is 
responsible for the ejecta north of the companion, drawn from the 
primary disk, and visible in the S 4 G image. During the prolonged 
interaction the tidal wave induced by the direct perturbation is 
able to propagate to the cen tral regions, consisten t with the inner 
spiral structure described by Zaritsky et al. (1993). 

has been very little testing of this theo r y with obser- 
vations. Modal studies by IBertin et al.l (1989a. b|) pre- 
dict that there should be amplitude variations along 
the arms from resonances between inward and outward 
moving spiral waves an d the precise amplitude varia- 
tions have been modele d (Elmcg reen fc Elmegreenl | 1984 : 
lElmegreen et all 119891 : lElmegreen fc Thomassonl 11993 : 
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Fig. 13. — Comparison of the m=2 Fourier amplitude profile of 
M 51 extracted from S 4 G 3.6fim image with the N-body simulations 
of Salo & Laurikaincn (2000a b). The simulation corresponds to 
the final snapshot of Fig. 1121 The outermost observed minimum 
at r Ri 180 arcsec is related to the change of the curvature of 
arms: in simulations a similar feature marks the region where the 
self-gravity sustained spiral wave separates from purely tidal outer 
structure. In the model the minima at r fs 40 and 20 arcsec are 
amplitude modulations, which arise due to interference between 
separate inward-propagating tidal wave packets. See also Figure 
11 for M51 at i? 2 5 « 300" 



IRegan & Elme greenl 119971) . Most recently, a new the- 
ory has been proposed for the formation of rings and 
spiral arms which argues in favor of chaotic orbits, con- 
fined by invariant manifolds ema nating from the LI fc 
L2 L agrangian points of a bar (|Romero-G6mez et al.l 
[2001 . Several comparisons between the results of this 
theor y and observations have b een carried out success- 
fully (|Athanassoula et al"]l2009bl lal) and several more will 
be possible with the S 4 G database. S 4 G will cover the full 
range of galaxies with varying spiral arm strengths and 
galactic structures. In Figure ITT] we show a preliminary 
result of the type of analysis that is possible with the S 4 G 
data. In this figure we show the arm-interarm contrast 
measured from the 3.6/im S 4 G images for NGC 7793, 
a flocculent galaxy, in contrast with M51, a grand de- 
sign spira l galaxy (see a model for amplitude variation 
of M51 in Salo & Laurikaincn 2000a). The contrasts are 
measured beyond R25, further than in previous optical 
studies. The spiral structure traced in the 3.6 /zm band 
is different in detail from the optical images, but still 
shows weak arms in NGC 7793 and increasingly stronger 
arms in M51. We will carry out a detailed analysis of 
spiral arm amplitudes and variations and compare the 
observational data to the models. 

For the NGC 5194 (M51) system, there may well be 
a connection between the spiral structure and the pres- 
ence of the compa nion (NGC 5195) . Indeed, detailed 
HI observations bv lRots et al.l ([1990D show the presence 
of a long tail, which emanates from the west side of 
the galaxy, running south of the spiral arm visible in 
the S 4 G image, and turning up towards the north fur- 
ther eastwards. The HI kinematics are very difficult to 
model with an inter a ction in volving a single passage, but 
iSalo fc Laurikainenl (|2000al Tbl) propose a model with a 
second passage (Figure fT2|). which gives a better fit to 
some of the aspects of the observations. The inner spi- 
ral structure in the main spiral is then amplified by the 
action of the companion, and wave interference causes 
amplitude modulations, which are well reproduced in the 



S 4 G observations (Figure IT5|) . This shows the interest of 
the S 4 G data for detailed modeling of the dynamics of 
individual galaxies. 

6.6. Early Type Galaxies & Dwarfs 

S 4 G contains 188 galaxies of T-type < and there- 
fore constitutes one of the largest and most homogeneous 
mid-infrared data sets for early type galaxies. The shapes 
of ellipticals and lenticulars may reflect different evolu- 
tionary paths, with wet mergers leading to disky shapes 
and dry mergers leading t o boxy profiles (iNaab et al.l 
[200l iPasauali et al.l [20071: iKang et alj 120071 ). In addi- 
tion, iKormendv et al.1 '{ 20091 ) shows that deep imaging 
can reveal departures from the Sersic profiles in ellip- 
tical galaxies, which are diagnostics of their formation. 
They found a dichotomy in which ellipticals that have 
cuspy cores at their inner radii can be separated from 
those which show an excess of light at the center. The 
S 4 G data offer us a unique opportunity to revisit the fre- 
quency of the different shapes and profiles as a function 
of stellar mass and environment, free of dust obscuration 
effects. 

Although dwarfs are the most abundant type of galaxy 
in the Universe, their low surface brightnesses of typically 
[iB > 22 mag arcsec -2 have restricted detailed observa- 
tions of these stellar systems. In particular, the charac- 
terization of the underlying stellar structure in these faint 
galaxies still remains largely unexplored. Using an abso- 
lute magnitude criterion of Mb > —17, we have identified 
a preliminary sample of 206 dwarf galaxies in S 4 G . The 
sample includes 87 late-type DDO dwarfs previously un- 
observed with Spitzer, to complement the 48 which have 
been observed. The S 4 G dwarfs therefore constitute a 
representative sample of nearby dwarf galaxies in which 
to study the properties of the underlying stellar compo- 
nent. 

Dwarf galaxies in S 4 G allow for a detailed analysis of 
the underlying stellar disk in both the 3.6 and 4.5/im 
data. By deriving the radial scale lengths for these galax- 
ies, more insight can be gained into the evolutionary link 
between dwarfs and their giant counterparts. At the 
same time, the relatively large sample of dwarfs in S 4 G 
gives a unique opportunity to look at how these param- 
eters may vary with environment. We can also derive 
stellar masses for the dwarfs which will be used to assess 
their stellar mass to dark-matter ratio and will likely 
provide crucial evidence to the longstanding debate of 

ly dark 

systems (see e.g., IStrigari et aI1[2008l) 

In addition to these areas, there are at least three other 
main areas of astrophysics that these data will address 
uniquely: the reported absence of bulges in some disk 
galaxies, the vertical stellar structure, and the Tufly- 
Fisher relationship. There is also likely to be significant 
amounts of spin-off science for S 4 G with its large and ac- 
curate inventory of mass and galactic structure in the 
nearby Universe, e.g., relationship between the Tully- 
Fisher relat ionship and propert ies of the fundamental 
plane (e.g., iZaritskv et al.l [2008). dust heating in outer 
disks, the interaction / merger fraction from a study of 
the debris around these galaxies, and diagnosis of AGN 
activity from tracing hot dust very close to the AGN 
from the [3.6]-[4.5] color. 
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7. DATA RELEASE AND TEAM POLICY 

The S 4 G team is committed to releasing all of the re- 
duced data upon publication of the entire S 4 G datasets 
from each of the pipelines. We expect the science-ready 
images from pipeline 1, the masks from pipeline 2, the 
galaxy properties and profiles from pipeline 3 and the 
input and output files from Pipeline 4 (including the 
mask, sigma images, images without NaN values and 
the GALFIT input and output parameter files) to be 
released on a staggered schedule using NED or IRSA for 
long term use by the community. Currently the prod- 
ucts of each of the pipeline are expected to be released 
in the supplement issue of the journals upon verifica- 
tion of their quality which will be on-going and occur 
as part of the various scientific investigations conducted 
by the team. We expect that some of the data and 
pipeline may need to be refined or enhanced in the fu- 
ture. As the data are still being acquired (data acquisi- 
tion is expected to be complete in 2011) and analyzed, 
our earliest release date for the the science-ready images 
(PI products) for the entire S 4 G data is not expected 
at least until 2012, and the supplement papers for the 
products from the other pipelines are expected to follow 
within a year after the Pipeline 1 supplement paper. In 
the meanwhile, to allow for maximum use of the data 
for science by the larger community, the collaboration 
has agreed on an open-door policy whereby any mem- 
ber of the astronomical community may join the team 
temporarily as a guest and pursue any type of science. 
The use of the data, publication policies and authorship 
on papers are clearly defin ed in our policy statement 
available here: http://www.cv.nrao.edu/^ksheth/ 
s4g/S4G_policy_v6.pdf . This statement is adopted by 
the S 4 G team and all of its guests. Interested members 
of the astronomical community are invited to read the 
policy and contact the PI or an S 4 G member to discuss 
the possibilities of joining the team. 

The noteworthy point of this policy is that our team 
has agreed to avoid carving out large science categories or 
areas to encourage the maximum use of the data within 
the team, and within the larger community, for any sci- 
ence that may be possible with these data. Our philos- 
ophy is that if multiple people are interested in a topic 
but want to work in separate teams then we will encour- 
age both teams to collaborate and address the scientific 
problems in independent but parallel investigations. 

8. SUMMARY 

S 4 G will provide the deepest, largest, and most ho- 
mogenous data set of nearby galaxies at 3.6 and 4.5/im. 
The spatial resolution will be unmatched at these wave- 
lengths, ~2", which is ideal for comparisons to ground 
based optical and near-IR observations, as well as 
MeerKAT, eVLA, ALMA and other future large radio ar- 
ray surveys. The sample of 2,331 galaxies includes most 
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morphological types and masses. The images will give 
an unparalleled view of stellar mass distributions and 
faint peripheral structures. All of the galaxies, as well 
as many additional galaxies from the Spitzer archives, 
will be reduced in the same way, giving data products 
that span the full range from science-ready images to 
measurements of global galaxy properties and individual 
components. S 4 G is designed to be a useful reference for 
many years to come. It will contain the best available 
mid-infrared data for individual galaxies and be among 
the most complete surveys for statistical studies. It will 
be useful to study the origin and evolution of galaxies 
and their dynamical components, and to supplement ob- 
servations of nearby galaxies at other wavelengths. 
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TABLE 1 

The S 4 G Sample. 



Galaxy 


RA (J2000) 


DEC (J2000) 


T 


M B 


Log(D 25 ) 


^Bcorr 




0.00784 




8.4 


-15.97 


1.06 


15.14 


PGC000143 


0.03292 


-15.46140 


9.9 


-15.83 


2.04 


10.00 


ESO012-014 


0.04511 


-80.34825 


9.0 


-18.17 


1.18 


13.72 


InCjC (814 


0.05418 


Id 1 /IKK /I 

16. 14004 


2.0 


-20.04 


1.64 


10.96 


UGC00017 


0.06198 


15.21822 


9.0 


-16.02 


1.39 


14.60 


NGC7817 


0.06635 


20.75182 


4.1 


-21.11 


1.52 


11.56 


ESO409-015 


0.09223 


-28.09915 


6.2 


-15.16 


1.02 


14.43 


ESO293-034 


0.10559 


-41.49592 


6.0 


-18.68 


1.40 


12.67 


NGC0007 


0.13916 


-29.91678 


4.9 


-18.22 


1.39 


13.18 


NGC0014 


0.14618 


15.81659 


9.8 


-18.25 


1.18 


12.35 


IC1532 


0.16464 


-64.37161 


4.0 


-17.89 


1.26 


13.86 


NGC7798 


23.99042 


20.74986 


4.1 


-20.20 


1.20 


12.56 


NGC7800 


23.99354 


14.80723 


9.7 


-19.45 


1.24 


12.61 



Note. — The full table is available in the electronic version of the 
journal. 



TABLE 2 
Sample RC3-like parameters 



Galaxy 






3.6/^m 






4.5/im 






sma 


cllip. 


PA 


sma 


cllip. 


PA 


NGC 0337 


25.5 
26.5 


90.6±1.5 
110.4±4.7 


0.287±0.018 
0.270±0.027 


118.4±2.0 
119.8±3.2 


84.1±1.3 
106.9±1.3 


0.278±0.037 
0.287±0.055 


118.5±4.5 
123.8±7.1 


NGC 4450 


25.5 
26.5 


209.0±6 
253.0±6 


0.293±0.016 
0.336±0.048 


1.4±1.8 
19.8±5.1 


199.9±... 
266. 7±... 


0.293±0.020 
0.418±0.026 


3.2±2.3 
18.1±2.3 


NGC 4579 


25.5 
26.5 


256.0±2.6 
278.6±17 


0.295±0.011 
0.251±0.018 


95.2±1.2 
78.6±2.4 


229.0±2.6 
263.4±12 


0.268±0.018 
0.282±0.027 


95.5±2.2 
98.3±3.4 



Note. — RC3-like parameters for three sample galaxies. For each galaxy and band 
we quote the semi-major axis (in arcseconds), the ellipticity and the position angle (in 
degrees) at two levels of surface brightness: 25.5 and 26.5 (AB) magarcsec -2 . The error 
in the semi-major axis is measured by noting the change in the corresponding radius 
when the surface brightness profiles are shifted up and down by the local measurement 
error plus 10% additional error, which is the maximum expected from any scattered light 
on the detector. For NGC 4450, the error cannot be measured due to contamination 
from a background star and is therefore set to the maximum error, 6", as discussed 
in the text. The errors in the ellipticity and position angle are from the ellipse fitting 
and reflect the rms in the formal fit of the corresponding isophotc. These errors are 
significantly larger than any error introduced by the shifting the profiles to measure 
the error in semi-major axis. 



TABLE 3 

Magnitudes and Concentration Indices for S 4 G Galaxies 



3.6/^m 

Galaxy mag (AB) C31 C42 



NGC 0337 11.432±0.001 
NGC 4450 9.630±0.001 
NGC4579 9.083±0.001 



4.5/im 
mag (ABj CJI 



2.67 2.74 
4.58 4.17 
5.29 4.54 



11.888±0.002 
10.112±0.001 
9.564±0.001 



2.82 2.90 
4.56 4.16 
5.60 4.71 
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